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We show that the effective field equations for a recently formulated polynomial affine model of 
gravity, in the sector of a torsion-free connection, accept general Einstein manifolds—with or with¬ 
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obtained from a gravitational Yang-Mills theory known as Stephenson-Kilmister-Yang theory. Ad¬ 
ditionally, we find a generalization of a minimally coupled massless scalar field in General Relativity 
within a “minimally” coupled scalar field in this affine model. Finally, we present a brief analysis 
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we prove that a Birkhoff-like theorem is valid for the analyzed sector. 
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I. INTRODUCTION 

A hundred years ago, A. Einstein presented by the first 
time the correct field equations of its relativistic theory 
of gravitation, known as General Relativity [I]. It was 
not only consistent with the available observations, but 
it also predicted several effects that have been consis¬ 
tently observed and measured—see Ref. [2] for a recent 
review—, being the last one the confirmation of the exis¬ 
tence of gravitational waves [3]. Thus, General Relativity 
has become the most successful theory of gravitational 
interactions nowadays. 

However, General Relativity seems to be an effective 
theory of a, yet to be discovered, more fundamental de¬ 
scription of gravitational effects. A first hint was pointed 
out by E. Cartan, after remarking that the curvature is 
not built up from a metric field but from a connection, 
and that a part of it was left behind if one chooses to work 
with the Levi-Civita connection [4-7] . One could also see 
this by considering Palatini’s approach in General Rela¬ 
tivity, in which both metric and connection fields are in¬ 
dependent, and noticing that generically the theory pos¬ 
sesses terms which do not appear in General Relativity 
(when it is coupled to matter). In addition, we mention 
a few aspects related with the quantization: Eirst, the 
fact that within the standard formalism of quantization. 
General Relativity is non-renormalizable [8-11], suggests 
the theory can be thought as an effective theory. Sec¬ 
ond, the Wheeler-DeWitt equations are not well-defined 
in general due to their nonpolynomial dependence on 
the Arnowitt-Deser-Misner (or ADM) variables [12-17]. 
One of the most successful quantization methods of grav¬ 
itational interactions—developed from ideas by Gambini 
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and Trias [18, 19], by Ashtekar [20, 21]— , hints that a 
first order formulation of General Relativity is not as 
complete as expected, since the Lagrangian formulation 
should include an extra term known as Holst term [22]. 
Finally, it is well known that once matter couples to grav¬ 
ity (in particular fermionic matter), the assumption of 
vanishing torsion in standard General Relativity is no 
longer preserved as a field equation, and the best ap¬ 
proach is to consider Cartan’s generalization of General 
Relativity [23-26]. 

Standing on the argument of the necessity of a quan¬ 
tum theory of gravity, several generalizations of Ein¬ 
stein’s General Relativity have been proposed. One im¬ 
portant class of these generalizations considers the con¬ 
nection field to be independent of the metric, known as 
Palatini type theories. Within this class of theories one 
finds the generalization due to Gartan—which considers 
a connection with nonvanishing torsion [4-7]— , but also 
a broader kind of generalization which do not require the 
metricity condition, see for example Refs. [24, 27, 28]. 

Interestingly, gravitational theories with a generic lin¬ 
ear connections have degrees of freedom which do not 
correspond to the expected massless graviton [29] . Albeit 
these gravitational theories with an arbitrary connec¬ 
tion possess ghosts, there are examples of well-behaved 
(stable) systems which have ghosts, such are the cases 
of the Pais-Uhlenbeck oscillator [30-33], or the higher 
derivative supersymmetric quantum mechanical system 
presented in Ref. [34]. The last sentence, allows to argue 
that a model cannot be disregarded for the presence of 
ghosts model, without a rigorous analysis of its stability. 

In order to undertake the premises of incompleteness 
and non-polynomial structure of the theory, in a pre¬ 
vious paper [35], we presented a polynomial (purely) 
affine gravitational model in four dimensions built up 
entirely on the basis of full diffeomorphism invariance. 
Worth to mention is that the most general action built 
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this way contains only power-counting renormalizable 
terms (which by no means guarantees the renormalizabil- 
ity of the system). We showed that its non-relativistic 
limit, around a homogeneous and isotropic spacetime, 
yields a geodesic deviation corresponding to the gener¬ 
ated by a Newtonian gravitational potential, despite the 
existence of torsion. We also highlighted the possibil¬ 
ity of using standard methods to quantize the model, 
due to its polynomial structure (unlike the earliest pro¬ 
posals [36, 37]), and also the likelihood of avoiding the 
uniqueness of diffeomorphism-invariant states from Loop 
Quantum Gravity programme [38] , due to the absence of 
a fundamental metric, or in other words, the absence of 
flux operators. 

Interestingly, by construction, the polynomial affine 
gravity proposed in Ref. [35] has no explicit terms leading 
to three-point graviton vertices, since all graviton self¬ 
interaction is mediated by non-Riemannian parts of the 
connection. This feature might allow to bypass the gen¬ 
eral postulates supporting the no-go theorems stated in 
Refs. [39, 40], where it was proven that generic three- 
point graviton interactions are highly constrained by 
causality and analyticity of the S'-matrix, and the only 
acceptable structure of the three-point graviton vertices 
is the one coming from General Relativity. 

The aim of this work is to show that in the vanish¬ 
ing torsion sector, the field equations of the polynomial 
affine gravity are a generalization of the Einstein’s equa¬ 
tions. The paper is organized as follows. In Sec. II, we 
present briefly the polynomial affine gravity model, and 
summarize the results obtained in Ref. [35] concerning 
the Newtonian limit of the model. In Sec. Ill, we restrict 
ourselves to the sector of the theory with vanishing tor¬ 
sion, and find the field equations; which surprisingly are 
a known generalization of the Einstein’s field equation. 
Then, in Sec. IV we show that under certain considera¬ 
tions the coupling of the gravitational model with a scalar 
field represents a generalization of the minimally coupled 
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system Einstein-Klein-Gordon. In Sec. V, we present 
a brief analysis of the propagators of the model, which 
allows us to count the propagating degrees of freedom. 
In addition, we corroborate that result with a method 
equivalent to the one of Dirac [41]. Finally, in Sec. VI 
we summarize and discuss the possible implications of 
the results. In order to make the paper self-contained, 
we have included three appendices: In Appendix A, we 
describe a dimensional analysis strategy that allowed us 
to build up the most general action consistent with dif- 
feomorphisms invariance. In Appendix B, we show a col¬ 
lection of known solutions to the effective field equations 
of the model. Finally, in Appendix C, we demonstrate— 
restricting ourselves to the vanishing torsion sector—the 
equivalent of the Birkhoff’s theorem of the theory. 

II. POLYNOMIAL AFFINE GRAVITY 

First of all, we highlight that the model constructed 
below has as only fundamental field an affine connection, 
and no metric field is needed neither for contracting nor 
lowering or raising indices. Moreover, in order to guaran¬ 
tee the correct transformation of the Lagrangian density, 
the geometrical objects used to write down the action 
will be the curvature and torsion of an affine connection, 
r^pcr, which accepts a decomposition on irreducible com¬ 
ponents as 

= f^(p.)+f'^[pa] = (i) 

where T^po- = r^(po-) is symmetric in the lower in¬ 
dices, Ap is a vector field corresponding to the trace 
of torsion, and jg ^ Gurtright field [42], satisfying 

and = 0. Since no metric 

is present, the epsilon symbols are not related by rais¬ 
ing (lowering) their indices, but instead we demand that 


Using the above decomposition, the most general action preserving diffeomorphisms is 


S'[r,r,A] = J d^a: 




+ S5 Rp,PpT-’P’'A, + Cl RppP.VaT’^'P^ + C2 Rp.PpVaT'^’P'' + D, 

+ D2 + S3 + S4 ixA,)ep,pa 

+ Ds T^'P’'V[xT^’P'^A,^ep,pa + Sg A,Vf^xAp) + Dy AxV[pA,] 

+ Si WipTP^P'^W^)T^A^ep,x. + E2 + T^AlT^,r,.rj.X,puj.p,ar(^p^ eis-,r,.ea.pp.esxar 

+ F2 epxp.e^ppueo^sar) + S3 ArCaP^xe^upa + S4 


( 2 ) 


In Eq. (2), we have dropped all terms which can be re- topological invariant (e.g. the Euler density). A straight- 
lated through partial integration, and those which are 
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forward dimensional analysis shows the impossibility of 
other contributions into the action, see Appendix A (in 
particular Sec. A 2). Interestingly, the above action turns 
out to be power-counting renormalizable, which does not 
guarantee renormalizability, but is a nice feature. The 
structure of the model yields no three-point graviton ver¬ 
tices, which might allow to overcome the no-go theorems 
found in Refs. [39, 40]. It is worth mentioning that all 
coupling constants are dimensionless, which might be a 
hint of conformal invariance of the model [43] , and further 
analysis of this invariance will be presented in Ref. [44]. 

Although there is no obvious equivalence between the 
action in Eq. (2) and General Relativity, particularly 
due to the lack of a source metric field, we analyzed in 
Ref. [35] the scalar perturbations and studied perturba¬ 
tive inhomogeneous sources to the connection field equa¬ 
tions, by considering a generic matter action. In the rest 
of this section we summarize the procedure which leads 
to the Newtonian limit of the model. 

First, we considered a static, homogeneous, and 
isotropic expansion of fields. 


— S^A a 


Ai) 


= SXoT + 


( 3 ) 




Albeit the inclusion of matter in a non-metric space- 
time is a nonstandard procedure, we assumed that the 
matter action would depend on the most general, sym¬ 
metric (p)-tensor density, gP’', built up with the available 
fields. The dimensional analysis in Sec. A 1 shows that 
this tensor density contains the same terms as what we 
call the “Eddington’s metric density” [45] , defined as the 
variation of the action with respect to the symmetric part 
of the Ricci tensor (see Refs. [36, 37, 46]), i.e., 

^—s[r]^^gp'^=rr (7) 

on{nv) 


From the geodesic equations and the linearized action 
on scalar perturbation, we found that only a few terms in 
the scalar perturbation affect the geodesic equation [47]. 
The effective geodesic deviation is given by 


T^i ^ '9'^Matter 

i 00 — — -- V 


Stt dg 


00 



( 8 ) 


is the usual Newtonian force induced by a massive source, 
where we restricted ourselves to consider gP'^ = gp''. 


III. RELATIVISTIC LIMIT 


with 5^^ = Next, by analysing the first 

order perturbations of the actions, we found a choice of 
the coupling constants which allows nontrivial solutions 
of these field equations, together with the non-relativistic 
limit of the geodesic equation. Then, we substitute the 
field components by their scalar perturbation decompo¬ 
sition, 

F (4) 

tP^^P ^ 

( 5 ) 

and 



+ (yiS^n + ( 6 ) 

+ 

+ (e™«n5p + + d^dr,dj,Zi}j , 

where the “L” sub-index identifies the longitudinal de¬ 
grees of freedom. 


Despite we obtained the correct Newtonian gravita¬ 
tional potential in the non-relativistic limit, see Eq. (8), 
a post-Newtonian approximation is necessary to explain 
gravitational phenomena like Mercury’s perihelion and 
the bending of light by a gravitational source. There¬ 
fore, we should go beyond the Newtonian limit in order 
to compare our model with General Relativity. 

Even though a comparison between the models does 
not require the imposition of vanishing torsion, for the 
sake of simplicity, at this stage we shall focus on a sector 
of the theory in which the connection is torsion-free. No¬ 
tice that the vanishing torsion limit—equivalent to take 
T^^p^ _>• 0 and A^ —>■ 0—cannot be taken at the action 
level, but in the field equations. Interestingly, the field 
equations from Eq. (2) can be consistently truncated un¬ 
der such requirements. It can be realized as follows, since 
only the terms of the action linear in these fields will be 
relevant, one can consider the effective action linear in 
the torsion’s fields, i.e., 

S-eff = J d^x (Ci Ra/.Vp + C2 Ryx^.)T''’PP, (9) 

and the only nontrivial field equation after the limit will 
be the one for the Curtright, T'^’Pp, 

= 0 , ( 10 ) 

with K a constant related with the original couplings of 
the model. We assume that the connection is compati¬ 
ble with volume form, i.e., it is equi-afhne [48, 49]. The 
equi-affine condition assures that the Ricci tensor of the 
connection is symmetric, and the contraction of the last 
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indices vanishes, thus the second term in the field equa¬ 
tion is absent and the gravitational equations are 

^[pR^l\u = ^■ ( 11 ) 

Equation (11) is a generalization of a condition known 
in Riemannian geometry as covariantly constant Ricci 
curvature—aka parallel Ricci curvature—, VpR^v = 0. 
All Einstein manifolds, whose Ricci tensor is proportional 
to the metric, Rpi, oc satisfy the parallel Ricci condi¬ 
tion due to the metricity condition, and therefore every 
vacuum solution to the Einstein’s equations solves the 
(simplified) field equations of our model. Consequently, 
the fact that the non-relativistic limit of the gravitational 
potential in Eq. ( 8 ) yields a Newtonian potential seems 
clearer, and we can argue that even the post-Newtonian 
corrections coming from General Relativity are present 
in the chosen scenario of our model. An additional com¬ 
ment, the parallel Ricci curvature assures that although 
the manifold is not Einstein in general, the metric has to 
be locally a product of Einstein metrics [50] . 

The Eq. (11) is related through the second Bianchi 
identity to the harmonic curvature condition [51], 

VaR^.% = 0, (12) 

and it can be shown that a manifold with harmonic curva¬ 
ture is equivalent (in four dimensions) to a manifold with 
harmonic Weyl tensor and constant scalar curvature [52] , 
or in other words the Ricci tensor is a Codazzi tensor [53] . 
For proofs of these equivalences, see Refs. [50, 54]. 

Notice that Eqs. (11) and (12) accept a geometrical 
interpretation equivalent to that of the field equations 
of a pure Yang-Mills theory, which in the language of 
differential forms are 

DF = Q, £)*F = 0, (13) 

where F = DA is the field strength 2-form (the curvature 
2 -form of the connection in the principal bundle, see for 
example Ref. [55, 56]), and the operator * denotes the 
Hodge star. Now, these Yang-Mills field equations are 
obtained from the variation of the action functional 

^YM = yTr(F*F), (14) 

and the Jacobi identity for the covariant derivative. 

Interestingly, the Eq. (11)— equivalently Eq. (12)— can 
be obtained from an effective gravitational Yang-Mills 
functional action [57-59], 

= J Tt{R*R)= J (15) 

where R £ M ®TM) is the curvature two- 

form, the operator * denotes the Hodge star, and the 
trace is taken on the bundle indices (see Ref. [55]). 

The gravitational model described by the action in 
Eq. (15) is called Stephenson-Kilmister-Yang (or SKY 


for short). However, the standard interpretation of this 
model requires a metric tensor. A quick analysis of the 
field equations a la Palatini shows that Eqs. (11) and (12) 
are obtained only from the variation with respect to the 
connection. Albeit this effective theory has been widely 
studied, the extra field equation (for the metric field) 
strongly constraints the solutions of the gravitational 
Yang-Mills [60]. However, that complication does not 
affect our affine interpretation of the theory, where the 
metric is absent. 

Intriguingly enough, the effective theory is not 
renormalizable—according to the arguments in Refs. [39, 
40]— , although it has been discussed in Ref. [61] that 
at least the ghost problem can be avoided through the 
arguments exposed in Refs. [31, 62-64]. 

IV. POLYNOMIAL AFFINE GRAVITY 
COUPLED WITH A SCALAR FIELD 

Until now, the most general diffeomorphism-invariant 
and power-counting renormalizable (gravitational) the¬ 
ory for an affine connection has been built, and we have 
showed that in a certain sector it is a generalization of 
General Relativity. In what follows, we show an at¬ 
tempt of including scalar matter into the model. The 
theory does not require a spacetime metric, but instead 
we should consider the most general, symmetric (p)- 
tensor density, built with the available fields (see 
Appendix Al), and use it to build Lagrangian densities 
for the matter content. Following our precept, the mat¬ 
ter content—scalar matter—should couple to given 
by 

gA*- = aVAT'^’"^ (16) 

with a, P and 7 arbitrary coefficients. 

We consider the action provided by the “kinetic term” 

S<t> = - J (a + /3 

which makes a nontrivial contribution to the field equa¬ 
tions once we restrict to the sector of interest. 

The equation for the Curtright field when the scalar 
field is turned on is (without lost of generality we fixed 
the coefficient Ci = I) 

- C^^i^Rpo^p - oS/[„{dp](i)du(l^ = 0 , 

which under our considerations simplifies to 

^[aRp\v-aV[a{dp\(l>dy(j^ =0. (18) 

In that case, we find a particular solution of Eq. (18), 
R^i/ 0id^(l)di/(l) — 
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where we have used a suggestive notation by denoting 
with a covariantly constant, invertible, and symmet¬ 
ric two-tensor. Notice that a is the gravitational cou¬ 
pling constant, which in Einstein-Hilbert gravity is pro¬ 
portional to Gn- This equation can be written in the 
more conventional form 

Rtj.u - ]^gfiuR + ]^g^,u{d(j)fY (19) 

Additionally, the second Bianchi identity imposes 

= 0. (20) 

This condition is, in the sense argued in Ref. [65], the 
equation of motion for the scalar field. Notice that, the 
Euler-Lagrange equation of motion for the scalar field 
yields no information after taking the vanishing torsion 
truncation. 

It is well-known that this system of equations can be 
obtained effectively from the Einstein-Hilbert action cou¬ 
pled minimally to a massless scalar held 

SeS=^J d^x (^ + 2A - . (21) 

We have show, so far, that our model of polynomial 
purely affine gravity is compatible with General Rela¬ 
tivity in the sector of vanishing torsion, and in addition 
the compatibility can be extended to include the interac¬ 
tion with a free scalar held. The analysis of this coupled 
system could be enhanced to include a potential for the 


scalar held, however, we are not considering this case at 
the moment. 


V. ANALYSIS OF PROPAGATORS 

A step toward the analysis of renormalizability of the 
model should be the computation of the free propagators 
of the model around a given background. The choice of 
background is necessary in order to simplify the com¬ 
putations. Our hrst task is to expand up to second or¬ 
der in the helds around a maximally symmetric Rieman- 
nian background. Maximally symmetric spaces possess 
remarkably simple curvature tensors [66], 

Ru.\ = K[g,pS^-g^X)^ ( 22 ) 

and the Einstein’s held equations in vacuum 

Rfi^=Agf_i„, (23) 

require that AT = ^. Therefore, in notation of Eqs. (3) 
we have that A = T = 0, and the contribution due to the 
torsion and non-metricity enters as a perturbative effect. 
Thus, most terms in the action (2) do not contribute 
to the propagator analysis. Additionally, the hrst order 
perturbation of the curvature tensor is 

Rfiv p = ~ T pi' nl pfi- (24) 

Under these assumptions, the relevant terms of the ac¬ 
tion (2) are those accompanied by the couplings Bi, Ci 
and Ei, and the second order perturbation of the action 
yields 


S'(^)[7,t,a] = j 


d^x Ci( - A7“ psgcjt^^'^ + V„7“ ps - V„7^ ps ) 

+ + E 2 (V„a^ 

- Vaap ) -b [Bi + i B2^Aeapjsgr,\G°‘^t^^^ ~ ^ BijAaagpjt^' 


(25) 


0:7 


Longitudinal and transverse projectors can be introduced 
in order to decompose the action and specify in a more 
standard fashion the dynamics of Eq. (25). This can 
be done since we have the metric compatible with 
the connection, and = g^'^V^. Additionally we will 

denote by A = the covariant Laplacian. Let us 

dehne Pl^i, = the projector into logitudinal 

vectors, and similarly P±^i, = — Pl^u is the projector 

into transverse vectors. Therefore, transverse indices are 
dehned such that = 0. 


Thus, the perturbations can be split into their longi¬ 
tudinal and transverse components as 

Op = a_L^-f VpOL, (26) 

(27) 

and 

7 pi/= 7-L pi/+ V 7sp;y + 2V(^7 -b 2V ^(p7vi/)^28) 
+ 2V(^Vjy)7v'''' -b 2V''V(^V,y)7L, 






6 


where all the indices are transverse except those of the 
derivative the subindex _L indicates that it corre¬ 
sponds to the fully transverse projection of the tensor in 
all of its indices, the subindex L indicates that this is the 
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The perturbed action is then 


longitudinal projection of the tensor, while the subindices 
s, a and v specify whether the tensor is symmetric, anti¬ 
symmetric or a vector representation. 


5 = J d^x (^ - (^ + AA + (^((Ci - 2 € 2 ) 1 ^^ + + 2C2)K^^ 

+ |(2Ci + C2)A2A - iciA3)7,„ + ((Cl - 2C2)AA - 2 C 2 ^^)ix^ + 8 + ( - 2 C 2 A 2 

+ {—E 2 E E 3 — — i?4)AA^ a_La^ + ts^^a ^gC'iAA7^ p + —C2AAaL — —CiA^A^l ~ + A)A7s^^^ ^ , 

(29) 


where A = V“Vq is the Laplacian operator. This action 
sets the dynamics of two transverse-traceless tensors (or 
specific combination of tensors), as well as two transverse 
vectors and two scalars, which are associated to eighteen 
degrees of freedom. 

We have checked, using the method introduced in 
Ref. [41], that the formal counting of degrees of freedom 
for the effective action in Eq. (9) yields eighteen degrees 
of freedom. We discover during the implementation of 
these program, that our model presents an extra symme¬ 
try associated to the metric independence property. 

In addition, when we restrict ourselves to perturba¬ 
tions around a flat (Minkowski) spacetime, the above 
equation simplifies greatly, 


S= 


-Ci(A7,^,)*AC 


/' 

+ ((Cl - 2C2)V‘^AV,7 vm - 2C2A7l,^ 

- 8Cie,;3.i.V“Atf" - 2C2Aai^) *AtO , 


and it is clearer that there are only eighteen degrees of 
freedom. 


VI. CONCLUSIONS 

In this paper we have shown that the polynomial ajfine 
gravity presented in Ref. [35] is the most general (poly¬ 
nomial) model for an affine theory of gravity, see the 
Appendix A 2, consistent with invariance under diffeo- 
morphisms. The action in Eq. (2), is parametrized by 
a nineteen-dimensional moduli space—after the bound¬ 
ary and topological terms have been dropped—, which 
in some sense made us wonder about a possible duality 
between this model and the one proposed in Ref. [28]. 


The field equations obtained from the action (2) are 
a set of very complex coupled partial differential equa¬ 
tions. However, the system is consistent with the limit 
of and going to zero. This consistency is 

what is usually called a consistent truncation of the the¬ 
ory. In this paper we have restricted ourselves to such 
a truncation of our model, which constrained us to the 
torsion-less sector of the moduli space. Additionally, 
this truncation leaves only one nontrivial field equation, 
Eq. (10). Moreover, by considering torsion-free equi- 
affine connections—those preserving a volume element— 
or a Levi-Civita connection—in the cases where the man¬ 
ifold is metric—, the remaining field equation simplifies 
even further to Eq. (II). 

The equation (II) is a generalization of the field equa¬ 
tions in General Relativity, in the sense that, any Ein¬ 
stein manifold (whose Ricci tensor is proportional to the 
metric) satisfy the parallel Ricci condition, 

VaR;.. = 0. (31) 

Although all Einstein manifolds satisfy the parallel Ricci 
condition, the contrary does not hold. Therefore, in Ap¬ 
pendix B we present a summary of known (metric) solu¬ 
tions to Eq. (II). In addition, in Appendix C we prove 
that a metric ansatz with spherical symmetry solves the 
field equation only if it is static, i.e., a Birkhoff-like the¬ 
orem for our model. However, the non-linearity of the 
equations does not assure the uniqueness of the solution. 
Nonetheless, certain equations become linear if the func¬ 
tions A and B in Eq. (C4) are equal. Therefore, the 
unique solution is the (Anti-)de-Sitter-Schwarzschild so¬ 
lution. 

We also pointed that our field equation (II) can be 
obtained from an effective action, known as Stephenson- 
Kilmister-Yang (SKY) theory [57-59], which is built up 
in the spirit of a Yang-Mills theory for gravity. Despite 
the SKY theory has been criticized by many people, in¬ 
cluding C. N. Yang himself [60], because its field equa¬ 
tions are higher order differential equations for the metric 
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field, in our model the fundamental gravitational field is 
the connection. Equation (11) is a second order partial 
differential equation for the connection. Although we 
have presented a series of metric solutions of the effec¬ 
tive field equations (restricted to the torsion-free sector 
of the theory), the correct way to proceed is to propose 
a connection ansatz compatible with the symmetries of 
the problem. Using this procedure, we expect to find a 
non-metric solution to the field equation [67]. 

We would like to highlight that, since all gravitational 
effects described by standard General Relativity are con¬ 
tained in our model, particularly in the sector of metric 
compatible, torsion-free connections, there is enough lib¬ 
erty to start considering new cosmological effects coming 
from both nonvanishing torsion and non-metric connec¬ 
tions with vanishing torsion, beside the results reported 
in Ref. [ 68 , 69]. 

As a hrst step toward the coupling of matter with our 
polynomial affine gravity, we showed that there is a con¬ 
struct which might play the role of an inverse metric den¬ 
sity. Notice that such a tensor density does not neces¬ 
sarily satisfy a non-degeneracy condition as required for 
a metric, but even though it is an interesting question to 
ask under what considerations does this construct satisfy 
the properties of a metric. Using this tensor density, we 
proposed a kinetic term for a scalar field (under diffeo- 
morphisms), and show that the effective field equations 
are a generalization to those of General Relativity, when 
we restrict ourselves to the torsion-free sector. Despite 
the fact that in this sector the field equation for the scalar 
field collapses, and cannot be obtained from the variation 
of the action, using the algorithm proposed in Ref. [65], 
the usual scalar field equation is recovered in virtue of 
the second Bianchi identity. 

Finally, we analyzed the second order perturbation of 
the action in Eq. (2), around a maximally symmetric 
spacetime. From this analysis it is possible to conclude 
that, at perturbative level, the number of degrees of free¬ 
dom is eighteen, coming from a pair of transverse, sym¬ 
metric, rank-two-tensors, and a pair of transverse vec¬ 
tors. Moreover, we applied the Lagrangian equivalent of 
Dirac’s formalism for constrained systems (see Ref. [41]), 
and obtain the same number of physical states—it is a 
nice feature, because there are known physical systems 
for which the perturbative analysis of degrees of freedom 
do not coincide with the exact one [70]. 
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Appendix A: Dimensional analysis method 

We built the model using six ingredients, a Curtright 
a vector (A^), the covariant derivative defined 
with the Levi-Civita connection (V^), both Levi-Civita 
tensors [e^ivXp and and the Riemannian curvature 

(Ruv^p)- Since the Riemannian curvature is defined as 
the commutator of the covariant derivative, it is not an 
independent field, so it will be left out of the analysis, 
and only five ingredients remain. 

Our interest is in general to build tensor densities. 
Therefore, we need to account for the number of free 
indices, and weight density of these quantities. Denote 
by N (< 1 >) the operator which count the number (and po¬ 
sition) of indices of the field $, being positive (negative) 
for upper (lower) indices. Thus, we have that 

A(r^’"^) = 3 N{Ap) = -1 V(V^) = -1 

N{ep,xp) = -4 = 4 

Then, for a general expression of the form 

the indices counting yield 

^(ya^byCg dg...e) ^ 

with 

n = 3a—b — c — Ad+ 4:6 = 3a — b—c + 4£, (A2) 

where we defined £ = e — d, for the sake of simplicity. 

In the same spirit, we define an operator which counts 
the weight density, 

^ ^ W{Ap) = 0 W{Vp) = 0 
W{ep,xp) = -1 W(e^"^^) = 1 . 

Thus, the weight {w) of the general expression above is 
given by 

w = a + £. (A3) 

1. Inverse metric density 

Now, we illustrate the usefulness of the dimensional 
analysis by building the most general symmetric (p)- 
tensor density, which we call the inverse metric density. 
This particular case fixes n = 2 and w = 1. 

Equation (A3) can solved by choosing either a = 1 or 
£ = 1, which imply that Eq. (A2) is restricted to 

3-6-c = 2, 


TABLE I. Possible terms contributing to the inverse density 
metric. 


a 

b 

c 

£ 

1 

1 

0 

0 

1 

0 

1 

0 

2 

0 

0 

-1 

0 

2 

0 

1 

0 

1 

1 

1 

0 

0 

2 

1 


TABLE 11 . Possible terms contributing to the Lagrangian den¬ 
sity. 


a 

b 

c 

£ 

1 

3 

0 

0 

1 

2 

1 

0 

1 

1 

2 

0 

1 

0 

3 

0 

2 

2 

0 

-1 

2 

1 

1 

-1 

2 

0 

2 

-1 

3 

1 

0 

-2 

3 

0 

1 

-2 

4 

0 

0 

-3 

0 

4 

0 

1 

0 

3 

1 

1 

0 

2 

2 

1 

0 

1 

3 

1 

0 

0 

4 

1 


or 

4 - 6 -c = 2 , 

respectively. The former, yields the terms 

and (A4) 

no other contraction of indices is allowed due to the sym¬ 
metry. The latter, yields another possibility, 

(A5) 

Finally, one can check that any other choice to solve 
Eq. (A3)— by allowing negative values of i —make im¬ 
possible to solve Eq. (A2). Ergo, there is no other term 
in a symmetric (g)-tensor density built up with these 
helds. Explicitly, this general tensor density is written in 
Eq. (16). The summary of this analysis is presented in 
Table I. 


2. Lagrangian density 

The work of building up the most general scalar den¬ 
sity with our ingredients is as simple as before, but the 
analysis is much longer. Thus, we will show the proce¬ 
dure with lesser details than before. 


First, we are interested in a scalar density, which sets 
n = 0 and m = 1 in Eqs. (A2) and (A3), i.e., 

3a — b — c + = 0, (A 6 ) 

a + e = l. (A7) 

A possible solution of Eq (A7) is a = 1 and £ = 0. 
Such choice allows four possible solutions of Eq. (A 6 ), 
but only three of them are nonvanishing. These generate 
the terms of Eq. (2) whose coefficients were called B^, 
Bi, i? 5 , Cl, C 2 , Dq, £>7 and i? 2 - 

If one choose £ = 1 for solving Eq. (A7), all possible 
solutions for Eq. (A 6 ) are either vanishing or yield topo¬ 
logical terms. 

The case with a = 2 and I = —\ gives three possibili¬ 
ties to solve Eq. (A 6 ). These choices give the terms the 
Lagrangian density whose coefficients are Bi, B 2 , II 4 , 
ZI 5 , El, and F 4 . 

For 0 = 3 and £ = —2 there are only two possible 
solutions of Eq. (A 6 ), and these yield the terms of the 
action with coefficients Di, D 2 , D 3 and £ 3 . And also, 
the choice a = 4 and £ = —3 solves the equations, and 
yield the terms in the action with coefficients Fi and £2 ■ 

Notice that for values o > 5, it is not possible to solve 
both Eqs. (A 6 ) and (A7) simultaneously. Therefore, we 
conclude that Eq. (2) is the most general action built up 
with these helds. A summary of the choices for building 
the Lagrangian density is presented in Table IT 


Appendix B: Known metric solutions to the parallel 
Ricci equations 

In this Appendix, we will present a brief compendium 
of known (Riemannian) metric solutions of Eq. (11) in 
four dimensions. The results presented in this Appendix 
is not original, and can be found in Refs. [50, 54, 71-73]. 

a. Einstein spaces. — The simplest example of a 
manifold with parallel Ricci is a Riemannian Einstein 
manifols. This example was explained on the main text. 

b. Riemannian products. — The Riemamian prod¬ 
uct of manifolds with harmonic curvature. 

c. Conformally fiat manifolds. — These are solu¬ 
tions in four dimensions if their scalar curvature is con¬ 
stant. 

d. Warped products. — Compact warped products 
of the form {S^ x M, dt'^ + f {t) ds"^(g)), where {M,g) is 
an Einstein manifold with positive scalar curvature, R > 
0 , and / is a positive function on satisfying the dif¬ 
ferential equation 

f-^R = cf, (Bl) 

for a negative constant c. 

Twisted warped products of the form (R x M)/1 have 
parallel Ricci, if the Z-action on the product metric in¬ 
volves an isometry of g. 
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Compact warped products of the form 
(Ml X M2, ■ {gi X (72)), where Mi are two two- 

dimensional manifolds, Mi has constant Gaufiian 
curvature, Ki < 0, and M 2 is an Einstein manifold with 
scalar curvature R 2 = —2Ki. The function / : Mi —R+ 
is a C°° solution of 


VV = cf for 0 0. 


Many examples and explicit conditions of Rieman- 
nian metrics with harmonic curvature, and conditions for 
them to have parallel Ricci, can be found in Ref. [74]. 

Although the above are examples of manifolds with 
parallel Ricci within the context of “Riemannian” geome¬ 
tries, examples of connnections with parallel Ricci can be 
found in symplectic geometries, see Ref. [75]. 

I 

The nontrivial field equation for the ansatz (Cl) are 




2r2/3 






-2 r^/—/-r/^U-/ + 


1 


dr 


J dp-' ) 2 

20^2 




- + 2 A + 8^^(0)61161) 


= 0 


Even though four out of the six nontrivial field equa¬ 
tions are independent, the system can be solved uniquely. 
First, the component along the (gvX) = (2, 0, 2) yields 
/ = T{t) + R(r). Then, one notices that the equation 
along {fJ,vX) = (1,0,1) contains a term depending solely 
on t, which implies that T vanishes, i.e., T is a constant. 
From here, the equation for (fii'X) = (2,1, 2) gets simple 
enough to be solved and yields, R{r) = 1 — T + f + 
or equivalently 


fir) = l + -+Pr\ (C3) 

r 

with a and /3 integration constants. Finally, it can 
be checked that albeit the equation along the compo¬ 
nents {fi,iy,X) = (1,0,0) was not used, it is satisfied by 
Eq. (C3). 


Appendix C: Explicit example: Schwarzschild 
ansatz, and BirkhofF theorem 

As pointed out, all vacuum solutions of Einstein’s 
equations—with and without cosmological constant— 
satisfy the parallel Ricci condition trivially [51]. 
Nonetheless, an important result in General Relativity 
is the stability of the Schwarzschild metric, known as the 
Birkhoff theorem [76-79], which implies that the gravi¬ 
tational collapse of a spherically symmetric astrophysical 
object cannot emit gravitational radiation. 

In this section we show that the simplest ansatz of a 
static and spherically symmetric metric, which solves the 
Eq. (11) is necessarily the (Anti-)de-Sitter-Schwarzschild 
metric. This should be taken as a first step toward a 
Birkhoff theorem, despite the fact that a complete proof 
is far from been achieved. 

Consider first the simplest metric ansatz, whose line 
element is 


ds^ = —f(t, r) + 


dr^ 




■dn^ 


(Cl) 


/ - 4 - 2 + 2 ( 3 - 2 r/ ) ^/ 


dt dtdr 


dr 


dr ■ 


df‘ 




(0)6ll6l) 


^ (p f3—f -2 p + 2r^—f 

2rp[ ^ dr^-' ^ ^ df’ 


(C2) 


The above shows that the simplest time-dependent 
generalization of Schwarzschild metric is not a solution 
for the field equations of the considered affine model. Un¬ 
less, the spacetime metric is static and exactly that for 
a Schwarzschild’s exterior solution, with or without cos¬ 
mological constant. 

However, if one tries a metric ansatz with two different 
factors, say 

dr^ 

ds^ = -A{t,r)dt^+A— + r^dn^, (C4) 
B(t,r) 

the field equations are so complicated, that although 
hints on the solubility were found—and yield to Eq. (C3) 
for both A and B —, the uniqueness of the solution has 
not been proven. The field equations for the metric (C4) 
are 
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Vry/,-f^7ylA — 


4^2^3^2 


+ 2r^B^A^^A + 2r^B‘^A^^^^B-Qr^BA^^B^^B - r'^B^A-^B^^A 


dr 


drS 


df^dr 


- r^B^A 


dr 

d^ 

dtdr^ 


dt 


+ 6 r^A^^B - rBA^ —B^ - 2 r^B^A—B + rB^A —A^ - 2 r^BA^ — B - 2 rB^A^ —B 


dr 


dr 


dt dtdr 

d 


dt dtdr' 
d^ 


d ^d 


dt dr 


d 


B-2 r^B^—B—A - r^BA—B - 2 rB^A —B —A 


dr 


d 


dt J dt 


d 


dr 


d 


dt^ 


r^B^A^^B + 2rB^A^^B -AB^A^ + -2r^B‘^A^B 1 -^A 

dr‘‘ dr dt dt‘‘ j dr 


d 


d 


Zr^B'^A^^B -Ar^B^A^A + ArB^A^ - ^ .lu ru 


dr 


dr 


52 


1 


dr^ 


4^254^3 


2 r‘^B^A^ 


d3 


"^B^A^^^B-^A - r'^B^A^B^A^ + 2r^B^A^-^B-^A + AB^A^-^B + 6r^A^-^B^ 


dtdr dr 


dt dr^ 


dt dr"^' 


dt 


dtdr^ 
c 

dt 


A 


d ^d „ ,92 




d 


- ( r^B^A-^B^A- 2 r'^B‘^^A^ + 2 r^B'^A^A-Ar^BA^B"^ + 'ir^B'^A^B 

or or or Or^^ dt Ot^ J dt 


+ ( r^B^A^^B - 2 r^B^A^A 

dr dr J dtdr 




2v42 


' d"^ (92 d d 

rA^^:—B + rBA—^A + rA—B—A 
dtdr dtdr dt dr 




+ 


1 


4ri?2^2 


2 r^B^A^^B + r^B^A^B-^A - r^B^-^A^ 
dr^ dr dr dr 


-2r^BA^B + r^B^^B^^A -A{B^- B^)A^ (sfjl + sin\e)5ll5i) 


= 0 , 


(C5) 


I- 

which are very complex. Nonetheless, the symmetries in the four-dimensional indices have been partitioned as 

some of the components of these equations seem to favour /i = (a, A), then the general form of the connection is 

the relation A = B as solution. In addition, we found a 
very interesting reference in which a proof of the Birkhoff 
theorem is given for an Einstein-Cartan theory coupled 
with the gravitational Yang-Mills [80] , which support the 
hand-waving symmetry argument above. 

Hence, further analysis requires a new strategy. In¬ 
stead of considering a metric, a connection ansatz is pro- 


X^xBxC^ 


+ SXSl:'gBcM<^{x) + S\Sf^St-^N,ix), 


where 7“bc and j^bc are the Levi-Civita connection for 
the metrics gab and gAB respectively, and the factors are 


otco-va (.U1101UC.L1115 cx iiicuiiv., a, uixxxi cxiiocxuz, lo jyfa( \ ^ _ 1 „ama 2 p{x) J AT = 9 d n 

posed based in the structure of a spherically symmetric 2^ 

metric of the form. 


ds^ = gab{x) dx“ dx^ A-e^A^) gABi^) di^^, (C6) 


From the connection ansatz (C7) the Ricci tensor is 
- 2V,niV„ - 2 NmN„) 


+ {BmN + gMN'^aM°‘), 


(C8) 


where x™ represent the temporal an radial directions, 

while ip^ represent angular coordinates. Notice that and the equations of motion are 

I 


y[xRp]u = 5\5^Ky^gu[mdi}n + 2ng^[^Ni^{x) - 2V[,(lV^]IV„) 

+ gMN (vNaM<^ + ^giaM- - 2 Ni - 2 Va [m-Ni) - 2 M‘^NaNi'). 


(C9) 


We have used that gAB{^) is the metric of a two-sphere, and its curvature tensors are well known. In addition. 
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since the metric gab{x) is two-dimensional too, its curva¬ 
ture tensors are related with the metric by the relation 

Rab^d = '^{sigM-5lgad). (CIO) 

As argue in Refs. [81, 82], the solution depends on the 
nature of the surfaces described by p = const., which in 
our parametrization translate to the nature of Na- We 
first consider the condition g°‘^NaNi, > 0, which allows to 
set p as a good radial coordinate. Therefore, the general 
metric ansatz gab is diagonal, and the radial direction 
would be p, i.e., 

gab = P) dta ® dtfc p) Apa 0 dpf,. (Cll) 

With this ansatz, the nontrivial components of Eq. (C9) 
are 


+ SYtp = 0, (C12) 

9p7^-47^^-16a-Vti = 0, (C13) 

= 0, (CM) 


5p(A,e20 + 8 - (^15) 

- Ap-'^e^PYpp + = 0 . 

From equation (C12) it follows that, 

- 8/3-2 = Gi(p). 


In addition, since 


„2p _ 2p 

cc — o 




= e 


2p 


4^-2 + jr 


(C17) 


it follows from Eqs. (C14) that, 

+ 2/3-2 =G2 (p). (C18) 


Equations (C16) and (CIS) implies that all temporal 
dependence on time within TZ, T and enters through 
an additive common function, say B{t,p). However, the 
consistency of these equations with the value of TZ calcu¬ 
lated from the metric ansatz, in Eq. (Cll), sets this ad¬ 
ditive function to zero, B(t,p). Therefore, neither TZ, T 
nor /3-2 depend on t. By extension—from the definition 
of T —, a could depend on time through a multiplica¬ 
tive function. Nonetheless, this multiplicative function 
can be eliminated of the problem by a redefinition of the 
time coordinate. 

At this stage, we have proved that any spherically sym¬ 
metric (metric) solution of the Eqs. (11) is necessarily 
static, which is one of the main statements of Birkhoff’s 
theorem. 

This result, simplifies the Eq. (C5) but the equations 
are still non-linear. However, one can get ride off the 
non-linearity by setting A{r) = B{r). In this linearized 
regime, according to Eq. (C3) the general (metric) solu¬ 
tion is the (Anti-)de-Sitter-Schwarzschild metric. 

Although we restricted ourselves to the linearized 
regime of the field equations, it is possible to go a bit 
further in the nonlinear regime. As a preliminary re¬ 
sult [67] , we have found that using the Frobenius methods 
with a Laurent series expansion, the leading order of the 
solutions is three-fold in the radial coordinate. These so¬ 
lutions correspond to the pairs given by (a, /3)—a (resp. 
/3) is the leading order of the series expansion for the 
function A{r) (resp. B{r ))—equals to (-1,-1), (2,-2) 
and (0,-1). 
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